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Abstract 

We systematically analyze the impact of current and foreseen cosmic 
ray antiproton measurements on the properties of Dark Matter (DM). 
We find that: 1) The current data from PAMELA impose constraints 
on annihilating and decaying DM which are similar to (or even slightly 
stronger than) the most stringent bounds coming from Fermi gamma 
rays, for hadronic channels and with fiducial choices for the astrophys- 
ical parameters. 2) The upcoming data from Ams-02 have the power 
to improve these constraints by slightly less than one order of mag- 
nitude and even to probe the thermal relic DM in the range 30—200 
GeV, for hadronic channels. However, with wider choices for the astro- 
physical parameters the uncertainty on the constraints spans between 
one and two orders of magnitude. We then explore the capabilities of 
early Ams-02 data to reconstruct the underlying DM properties in 
the case of a positive detection of a significant excess (attributed to 
DM annilations) over the background. For hadronic channels, we find 
that Ams-02 should be able to somewhat determine the DM mass 
and the cross-section, but not the specific annihilation channel nor 
the branching ratios. If other more exotic annihilation channels are 
allowed, the reconstruction will be more challenging. 
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1 Introduction 



Overwhelming evidence for the existence of Dark Matter (DM), in the form of an unknown 
particle filling the galactic halos, originates from many observations in astrophysics and cos- 
mology: its gravitational effects are apparent on galactic rotations, in galaxy clusters and in 
shaping the large scale structure of the Universe [1]. On the other hand, a non- gravitational 
manifestation of its presence is yet to be unveiled. One of the most promising techniques is the 
one of indirect detection, aimed at identifying excesses in cosmic ray fluxes (e ± ,j, u,p, d) which 
could possibly be produced by DM annihilations or decays in the Milky Way halo. The current 
experimental efforts mainly focus in the GeV to TeV energy range, which is also where signals 
from WIMPs (Weakly Interacting Massive Particles) are expected. 

Concerning e , the well known excesses detected by Pamela, Fermi and Hess have been 
indeed interpreted in terms of DM (see e.g. [2], among many) but then, given the tensions or 
plain inconsistencies with gamma ray and other constraints (see e.g. [3] for a concise review), a 
more mundane explanation in terms of some previously-unaccounted-for astrophysical process 
now seems preferred. In turn, this means that e ± lose their appeal as far as DM searches are 
concerned: the large excesses likely drown the possible signal from all but the most extreme 
(very large cross section, very large mass) DM models. 

In 7-rays searches, with a few exceptions, only constraints on DM are reported, from the 
observation of many different targets and in different bandwidths (see e.g. [4] for an up-to-date 
overview). One notable exception consists of the claim of detection of a 7-ray line at around 135 
GeV from the Galactic Center (and possibly elsewhere) in Fermi data. This feature, however, 
needs further scrutiny and we will not further address it here. 

In neutrino searches, only constraints are reported, both from observation of the Galactic 
Center or the Sun (see again [3] for references). 

Within this thriving area of research, we will focus on antiproton searches, for several 
reasons. First, antiproton production is relatively universal in DM models: as long as DM 
annihilates (or decays) into quarks or gauge bosons, p will be produced abundantly in the 
hadronization process, and even in leptophilic channels some p are generated, as we will discuss 
below. On the other hand, the yield of p from galactic astrophysical processes is relatively less 
important, a fact which makes for a potentially very good signal/background ratio. Second, 
the prediction of p flux from astrophysical processes suffers from a relatively small uncertainty 
(as we will recall later) and the error band is particularly narrow around tens of GeV, which 
is precisely a fortunate window for the emergence of a DM signal. Third, the propagation of p 
in the Galaxy is relatively under control and is generically less dependent on the details of the 
galactic geography then, say, e ± propagation. Finally, data on p, currently from Pamela, are 
exquisitely precise and even better is foreseen for Ams-02, so that the constraining power can 
be predicted to be very good. 

The aim of this paper is therefore to investigate the constraints and the projected sensitivity 
to DM signals in the p channel. A selection of past works with a similar aim, usually in some 
specific model or with the Pamela/Fermi/Hess e ± anomalies in mind, includes [5, 2, 6, 7, 
8, 9, 10, 11, 12, 13, 14]. A particularly thorough analysis has been done in [15] (with which 
we will compare later on). In this work we aim at a more systematical and model-independent 
analysis of antiproton capabilities than most previous works. A first step in this direction has 
been made in [16]: we here extend to several different channels and to a large range of DM 
masses. 
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The rest of this paper is organized as follows. In Sec. 2 we recall the basics of the production 
and propagation of antiprotons (from DM as well as from astrophysics) and the choices we make 
for the analysis. In Sec. 3 we present our results: we first compute the current constraints from 
Pamela data; we then derive the sensitivities expected from Ams-02 data (if no signal will be 
seen) and finally we investigate the reconstruction capabilities of Ams-02 itself assuming that 
a positive signal will be seen. Finally in Sec. 4 we conclude. 

2 Antiprotons from DM and astrophysics 

In this Section we briefly recall the basics of production of antiprotons in DM annihilations or 
decays ('primary antiprotons') and their propagation in the galactic halo of the Milky Way. 
We also specify our assumptions for the astrophysical p background ('secondary antiprotons'). 

2.1 Primaries 

The term 'primary antiprotons' refers to the p flux produced by DM annihilations or decays 
in the galactic halo. While we here briefly review the basic points focussing in particular on 
the choices we adopt, we refer for any other detail to the standard literature on the subject 
(e.g. [17] and references therein). 

We consider several different annihilation or decay channels, in a model independent way. 
We start by considering seven of them: 



In practice the spectral shape of the resulting p fluxes from all these channels are not very 
different, due to the qualitative fact that p are the final product of a complex hadronization 
process which 'washes out' possible distinguishing features. This fact, we anticipate, will imply 
that the channel cannot be easily determined on the basis of data fitting (see sec. 3.2.2). 
Moreover, some of the channels yield p spectra which are practically indistinguishable, both in 
terms of shape and normalization. This is the case for bb and it (of course for the range of DM 
masses in which both are kinematically allowed), for W + W~ and ZZ as well as for and 
t + t~). Therefore, for all practical purposes, only four of the channels lead to actually different 
results (bb,W + W~ , fi + fi~ and 77) and thus we will mainly focus on these in the following. 
Along these same lines of argument, the results for other channels that we do not consider can 
be safely deduced from channels with similar spectra: for instance, the hh channel (where h 
denotes the ~ 125 GeV SM higgs boson) is equivalent to the bb one, the e + e~ is equivalent to 
the one etcetera. 

Notice that, alongside with the 'traditional' quark and weak gauge boson channels (col- 
lectively: 'hadronic channels'), we consider leptonic channels (/z + /i~, r + r~) and a direct 77 
channel. Antiprotons are produced in these channels because we take into account electro- 
magnetic and electroweak corrections (i.e. splitting of the initial particles and radiation of 
electro(weak) bosons) to the tree level production process. Not surprisingly, the yield of p from 
these processes is suppressed and therefore (we anticipate) the resulting constraints will be 
much weaker than for hadronic channels. Nevertheless, we include these channels for complete- 
ness and for the interest that they have carried in recent years for the interpretation in terms 
of DM of the Pamela and Fermi anomalies and of the 135 GeV line in Fermi data. 
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Antiproton propagation parameters 
Model 6 ICq [kpc 2 /Myr] V conv [km/s] L [kpc] 
MIN 
MED 
MAX 

Table 1: Astrophy steal parameters. Left panel: the parameters for the considered DM profiles, to 
be plugged in the functional forms of eq. (2). These specific values are derived as discussed in [17]. 
Right panel: propagation parameters for anti-protons in the galactic halo (from [22, 23]). Here 5 
and fCo are the index and the normalization of the diffusion coefficient, V conv is the velocity of the 
convective wind and L is the thickness of the diffusive cylinder. 



The antiprotons produced by annihilations or decays in any given point of the galactic DM 
halo are then subject to propagation in the galactic environment. The prediction for the fluxes 
collected at Earth depends therefore on the assumed distribution of DM in the halo and on the 
assumed propagation parameters. Concerning the DM halo profiles we consider three of the 
usual possible choices motivated by numerical N-body simulations or by observations, namely: 
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The Einasto [18, 19] and NFW profiles [20] (favored by simulations) are peaked towards the 
Galactic Center while the Burkert profile [21] (possibly more in agreement with observations of 
galactic rotation curves) sports a constant-density core at the center. The specific parameters 
to be plugged in eq. (2) are given in Table 1. 

Concerning the propagation parameters, we adopt the standard three sets reported in Ta- 
ble 1 (right), derived from consistency with ordinary Cosmic Ray (CR) measurements. These 
parameters are to be plugged into the equation that governs the diffusion-loss propagation of 
antiprotons in the galactic environment: the set denoted by 'MIN' ('MAX') corresponds to a 
minimal (maximal) final yield of p and therefore they bracket the uncertainty. We will briefly 
comment later on how favored each set should be considered. 

Concretely, we obtain the p fluxes at Earth (post-propagation) from the numerical products 
provided in [17]. These fully take into account the EW corrections mentioned above [24]. As 
an example, the p flux from a model we will consider below is plotted in fig. 1, right. During 
the whole analysis, only kinetic energies above 10 GeV are considered, to minimize the effect 
of solar modulation. 



2.2 Secondaries 

The term 'secondary antiprotons' refers to the smooth p flux produced by standard CR pro- 
cesses (essentially the collisions of primary p on the interstellar medium), which constitutes the 
background to the DM signal. 

We do not calculate such background ourselves, but use instead the recent determination 
obtained by [15]. There, the two extreme cases are the Kolmogorov (KOL) and Kraichanian 
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Figure 1: Examples of antiproton fluxes and data. Left Panel: astrophysical p background, 
superimposed to the current data from Pamela and to one of our realizations of mock data from the 
Ams-02 experiment. Right Panel: the same but including a contribution from DM annihilation giving 
rise to an excess on top of the background which could be detected by Ams-02. The example chosen is 
DM annihilating into bb with todm = 1 TeV and (av) = 5 x 10~ 25 em 3 s for an Einasto halo profile 
and MED' propagation parameters. 



(KRA) models: we choose a spectrum which goes roughly through the middle of these (for 
illustration, it is the one plotted in fig. 1). We will refer to this in the following as the 'fixed 
background'. 

We then have to model the uncertainty of such background function: we let the amplitude 
vary by 10% and the spectral index by ±0.05. This envelops fairly generously the different mod- 
els and the astrophysical uncertainties. Technically, we just multiply the fixed background by 
AT P (where T is the antiproton kinetic energy), and allow A G [0.9, 1.1] and p e [—0.05, 0.05]. 
In the following, we will marginalize with respect to the A and p parameters, i.e. we will identify 
the values that give the best fit for a given assumed signal and a given set of data. We refer 
to this procedure as 'varying background' (see the detailed discussions below). In absence of a 
signal, the best fit background with respect to the Pamela data is obtained with an amplitude 
of A = 1.08 and a slope p = —0.05: it gives a Xo(AbR>) = 4.62 for 8 data points (we always 
restrict to kinetic energies T > 10 GeV). 

We note that, in previous analyses (e.g. [16]), the background used was the one in [25] with 
wider intervals for the variating amplitude (40%) and slope (±0.10). The obtained uncertainty 
band was thus broader and these choices would allow to fit the Pamela data slightly better. 
However, we decide to make use of the function deduced from [15] since this analysis takes 
recent cosmic ray and proton spectra into account and the uncertainties are clearly smaller. It 
is of course possible to get the same larger bandwidth from the new background by increasing 
the amplitude and the slope. 

Before moving on, we comment on the robustness of the assumption that the background 
follows a smooth shape, essentially consisting of a power law (above 10 GeV). While this is 
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Figure 2: Annihilating DM: current constraints. Left Panel: current constraints from the 
antiproton measurements by Pamela, for different annihilation channels. The areas above the curves 
are excluded. The dashed lines reproduce the j-ray constraints from [30], for the same channels. The 
symbols individuates the parameters used for the analyses in Sec. 3.2.2 while the horizontal band signals 
the thermal relic cross section. Right Panel: illustration of the impact of astrophysical uncertainties: 
the constraint for the bb channel spans the shaded band when varying the propagation parameters 
(dashed lines) or the halo profiles (solid lines). 



expected to follow from most astrophysical processes, it may not always be strictly the case. In 
fact, some astrophysical sources, such as supernova remnants [26] or plasma phenomena [27], 
could induce a different spectrum. However, such deviations from an effective power law are 
typically expected at energies higher than those of our interest [28]. Nevertheless, in case of 
a positive detection of a deviation, it will be crucial to keep in mind these possibilities when 
working on its interpretation. 



3 Results and discussion 

3.1 Current antiproton constraints from PAMELA 

The currently most precise measurement of the CR antiproton flux is provided by the Pamela 
satellite [29], as anticipated in the Introduction. The data, reproduced in fig. 1, extend from 
kinetic energies of less than 1 GeV to about 180 GeV (although we use only the portion above 
10 GeV to avoid the effects of solar modulation). 

The total antiproton flux is given by the sum of the DM and the astrophysical contributions, 

0totOiDM, (o-v); A,p) = 0dmOdm, M) + </>bk g (Ap)- (3) 

For fixed values of the DM particle mass m-£>M and the thermally averaged annihilation cross- 
section (<jv), the astrophysical background is optimized within the uncertainty bandwidth in 
order to minimize the \ 2 °f the total flux with respect to the data. During this procedure, 
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Decay constraints from antiproton flux 



Astrophysical uncertainties on the constraints 




Figure 3: Decaying DM: current constraints. Left Panel: current constraints from the antipro- 
ton measurements by Pamela, for different decay channels. The areas below the curves are excluded. 
The dashed lines reproduce the 7-ray constraints from [31]. Right Panel: illustration of the impact of 
astrophysical uncertainties: the constraint for the bb channel spans the shaded band when varying the 
propagation parameters (dashed lines) or the halo profiles (solid lines). 



the optimal values of the amplitude A opt G [0.9,1.1] and the slope p opt G [—0.05,0.05] are 
determined. Then, to find the 2a exclusion contour in the (moM, (&v ))-plane, the required 
condition is 

AxVdm, (av)) = X 2 (^dm, (o-v); A opt ,p opt ) - xl(A ,p ) < 4. (4) 

We show the resulting contours in figures 2 (for annihilating DM) and 3 (for decaying DM). 
In the left panels we have fixed the astrophysical parameters to a fiducial choice: the halo 
model is chosen to be Einasto and the propagation parameters to be 'MED'. The constraints 
on (av) reach the level of ~ 10 -25 cm 3 /sec for the quark or gauge boson channels, at a DM 
mass of about 100 GeV. For decaying DM, the constraints on the lifetime r reach ~ few x 10 27 
sec. The bounds for the exotic channels (77 and leptons) are much less stringent, as expected. 
The constraint from the leptonic channels is lifted below todm — 300 GeV, consistently with 
the fact that EW corrections are effective only for large 7Bdm [24]. The 77 constraint exhibits 
two regimes: at moM < 300 GeV some p are produced by 7 showering; at moM > 300 GeV 
EW corrections set in and increase the p yield, resulting in a constraint that follows the shape 
of the leptonic one. 

These results show that, for the hadronic channels, the constraints inferred from the data 
of Pamela are already competitive with (or even slightly better than) the bounds from Fermi 
7-rays (reported as dashed lines in figures 2 and 3, respectively from [30] and [31]), for both 
annihilating and decaying DM. This is a non-trivial conclusion, especially considering that the 
prospects of improvement in 7-rays are arguably less promising than in antiprotons (the Fermi 
satellite will at most double the statistics in the final years of its lifetime, while for antiprotons 
Ams-02 can in principle accumulate high precision data for several years to come). 
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On the other hand, the astrophysical uncertainties for antiprotons are non-negligible: the 
choice of propagation parameters and halo profile can change the constraints by one or two 
orders of magnitude, as illustrated in the right panels of fig.s 2 and 3. The impact of changing 
the propagation parameters from 'MIN' to 'MAX' is particularly important as it accounts by 
itself for more than a one order of magnitude shift in the bounds. 1 We also notice that, as 
expected, the impact of changing the halo profile is more important for the annihilation case 
than for the decay one. In fact, the differential p flux is sensitive to the halo density squared 
for annihilations while it is to the first power in the case of decay. 



3.2 Projected sensitivity of Ams-02 

In this section we move to study the reach of the upcoming early Ams-02 antiproton results 
on DM. 

We will therefore need to produce mock Ams-02 data and their associated error bars, which 
we do by implementing the following procedure. We adopt a linear approximation of the rigidity 
resolution of Ams-02 as given in [35]: 

AT 

r(T) = — = 0.0042 x T + 0.1, (5) 

where T is the kinetic energy of the incoming antiproton. This allows us to determine a realistic 
p energy binning. We assume that the p flux will be measured up to 300 GeV (we will use only 
the portion above 10 GeV, as already discussed, for a total of 27 simulated data points), thus 
extending the range of Pamela to higher energies. Next, we need to estimate the errors on the 
measurement of the antiproton flux. The number of collected and reconstructed antiprotons in 
Ams-02 in a given bin i centered around a kinetic energy Tj will be given by 

Ni = ea(Ti) 0^) AT t At, (6) 

where e is the efficiency, a the geometrical acceptance of the apparatus (which is in general a 
function of the energy), the antiproton flux, AT the width of the kinetic energy bin and At 
the exposure time. Assuming a Poisson distribution (AN = VN), the statistical error on the 
antiproton flux in bin % is therefore 



We assume systematic errors of 5% based on the analysis in [36] and we sum them in quadrature 
with the statistical ones. The p efficiencies of the different subdetectors of Ams-02 are expected 
to be close to 1 (or in any case well above 0.9) for the whole range of energies of our interest [37]. 
Therefore we set the overall efficiency e = 1. The modelization of the acceptance of the 
experiment is the most critical point and the one which is least accessible from outside the 
collaboration. We will therefore consider two options: 



1 It is worth noticing, however, that the propagation model 'MIN' has already been disfavored due to its 
small thickness of the diffusion cylinder. Studies on synchrotron emission, radio maps and low energy positron 
spectrum exclude a diffusion halo scale smaller than 2 kpc at 3cr level [32]. Moreover, a bayesian analysis on 
B/C, 10 Be/ 9 Be ratios and carbon and oxygen spectra suggests a thickness in the interval [3.2,8.6] kpc at 95% 
confidence level [33]. See also the analysis in [34], whose results also disfavors a I kpc thickness, although the 
observation region is not optimized and the results are weakened by the inclusion of Dark Matter. 
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1. At = 1 year and an acceptance a = 0.147 m 2 sr for 1 GeV < T < 11 GeV, but an energy 
dependent one (contained roughly in the interval a ~ 0.02 4- 0.04 m 2 sr) for energies 
T > 11 GeV, by interpolating the points in fig. 8 of [38]; 

2. At = 3 years with a constant acceptance for all energies a = 0.1 m 2 sr. 

The second case is obviously much more optimistic, as it gives rise to a much larger accumulated 
statistics. The first case might be more realistic though: in Ams-02, the identification of 
antiprotons at high energy will likely have to rely on the use of the calorimeter, thus reducing 
the pure geometrical acceptance of the experiment. We can expect the real data to lie in 
between these two scenarios. The resulting set of mock data (for option 1, for definiteness) is 
reproduced in fig. 1, both for a scenario with no excess (left) and for a scenario where an excess 
over background is detected (right). We indeed now move to discuss these two cases. 

3.2.1 Foreseen constraints 

First, we consider a scenario in which the measurements of Ams-02 do not shown an antiproton 
excess and the data points follow the background, as represented on Figure 1 on the left. 

We generate the mock data using the 'fixed background' spectrum function. We then 
compute the constraints for annihilation and decay of DM with the same method discussed 
in Sec. 3.1 (therefore, in particular, implementing the 'varying background' procedure). We 
believe this 'workflow' to be realistically similar to an actual data analysis. The results are 
presented in Figure 4. 

It is apparent that, with Ams-02, the constraints can be improved by about an order of 
magnitude or slightly less with respect to Pamela. For the bb channel (and for the tt one, when 
kinematically open) the bounds will reach the thermal value of the annihilation cross-section 
(2.7 x 10 -26 cm 3 /sec) in the range of masses 30 — 200 GeV. The more optimistic choice for the 
acceptance and the operating time (option 2 in 3.2) improves the projected reach by at most 
another half an order of magnitude or so, except at low masses (i.e. low energies) where the 
systematic errors dominate. 

We only show the curves corresponding to the 'Einasto, MED' choice of astrophysical pa- 
rameters: other choices would modify the bounds in the same way we discussed in Sec. 3.1. 
However, we also mention that a significant collateral improvement expected from Ams-02 will 
be to reduce the astrophysical uncertainties: using low energy spectra and the B/C ratio, as 
described above [32, 33], the propagation of cosmic rays should be better pinned down and the 
error band reduced. 

In [15] a similar analysis has been performed, for annihilating DM. Our results in fig. 4 
(left) compare almost directly with their fig. 17, albeit on a reduced range of masses and for 
the specific models that they consider. The constraints that we find are less stringent than 
theirs by about one order of magnitude. We attribute this to the fact that they assume a larger 
acceptance and, possibly, to their use of a larger range of energies. 

3.2.2 Reconstruction capabilities 

We now consider a scenario in which an identifiable excess is detected by Ams-02 in the antipro- 
ton flux and it is attributed to DM, and we study how well it will be possible to reconstruct 
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Annihilation sensitivities from antiproton flux Decay sensitivities from antiproton flux 




DM mass m DM [GeV] DM mass m DM [GeV] 

Figure 4: Annihilating and decaying Dark Matter: future sensitivities. Left panel: the 
sensitivity of early Ams-02 antiproton measurements for different channels. The solid lines assume 1 
year of data taking and a realistic acceptance, the dotted lines assume 3 years of data taking and an 
optimized acceptance. The rest of the notations are like in fig. 2. Right panel: the same for the case 
of decaying DM. 



the underlying DM properties. We focus only on annihiliating DM, for definiteness, in this 
Subsection. 

We need to assume values of the DM mass and annihilation cross-section which are still 
allowed by the Pamela data, but which would be probed by Ams-02. We therefore select 
three benchmark models: 

▲ m DM = 85 GeV (av) = 2.7 x 10~ 26 cm 3 s _1 , 
★ m DM = 300 GeV (av) = 1.5 x 10" 25 cm 3 s -1 , 

■ m DM = ITeV (av) = 5 x 10 -25 cm 3 s -1 , 

which are denoted by the corresponding symbols in Fig.s 2, 4 and 5. Each of these models 
represents a different region of interest in the (rriDM,(o~v))- i p\aiie: 

▲ The cross-section corresponds to the thermal annihilation cross-section, favored by cos- 
mological observations. Because of the relatively small mass, the DM signal affects mainly 
energies below 10 GeV (which are not considered for the analysis). 

-A- The p signal sits squarely in the energy range probed by Ams-02 . For the given mass, 
the cross-section is chosen at the limit of the exclusion contour of PAMELA. The analysis 
of this model should be the most straightforward. 

■ The DM signal starts to have an important contribution for energies around 100 GeV, 
where larger uncertainties are present. The lack of data for high energies should pose 
problem for the reconstruction of this model. 
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We assume an Einasto halo profile and we fix 'MED' propagation parameters. We sum the p flux 
resulting from these DM models to the 'fixed background' and we generate the corresponding 
mock data, plotted in Fig. 1, right. We thus determine a posteriori the regions of the parameter 
space which would be identified, at a given C.L., by a blind analysis of such data. 

Mass and cross-section for bb channel 

First we assume that the annihilation channel is fixed and known: XX ~^ bb with 100% Branch- 
ing Ratio (BR). In the left panels of Fig. 5 we show the regions identified at 90% and 99% C.L. 
(corresponding to A% 2 = 4.61,9.21 for 2 d.o.f.), with the ellipse corresponding to the 'fixed 
background' hypothesis and the extended shaded areas corresponding to the implementation 
of the 'variating background' procedure. 

-A- We focus first on the results for this model. We see that todm is determined fairly well, 
within 50% of the true value or so. For a fixed cross section, values of the mass which 
are smaller than the true value are more quickly disfavored than values which are larger 
than the true one. This is consistent with expectations since, for small masses, the signal 
will have its maximum at small energies, which are measured more precisely. 
We also notice that the ellipse corresponding to the 'fixed background' analysis displays 
the degeneracy between and (av) , since the product of these two quantities appears 
in the determination of the DM flux. On the other hand, the 'variating background' 
analysis destroys this simple analytical dependence. 

▲ Here, we display only the 90% C.L. contours: the 99% ones would artificially extend to 
very low DM masses, as a consequence of the 10 GeV cut imposed by solar modulation. 
From the point of view of the \ 2 i the points in the 99% contour would favor a pure 
background. The shape is unusual probably because of the cut at 10 GeV. 

■ The 99% C.L. region extends to very large masses as the signal moves into the region in 
which the measurements are least precise. In addition, as the signal affects large kinetic 
energies, any contribution with a mass > 1 TeV can fit the the data at 99% C.L. provided 
the cross-section is large enough. 

In all cases, the annihilating cross section is determined within an order of magnitude at best. 
Mass and the cross-section for different annihilation channels 

Next, we use still a single annihilation channel (xx ~^ bb with 100% BR) to produce the 
mock data and we investigate how the same signal could be interpreted in terms of a different 
channel. The results are shown in the central panels of Fig. 5, where the contours enclose 
the 90% C.L. regions for the indicated channels. Not surprisingly, interpreting the data in 
terms of annihilation into the W + W~ channel identifies a region which is very similar (a part 
for the kinematical cut in the ▲ case) to the one corresponding to the 'true' bb annihilation, 
since the p spectra and the yields of the two channels are very similar. On the other hand, 
interpreting the data in terms of annihilation into 77 (/i + /i~) would result in reconstructing a 
DM mass smaller (larger) than the true one and a cross section much larger than the true one. 
In Tables 3 and 4 (upper parts) we report the best fit values of DM mass and annihilation cross 
section (identified by a colored point in Fig. 5) for each channel. We also report, in the last 
column of the Tables, the A% 2 of each one of these points with respect to a pure background 
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Reconstruction of an Excess XX"* bb 



Reconstruction of an Excess at 90 C A C.L. 
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Figure 5: Reconstruction capabilities of DM properties by Ams-02. Left Panels: Reconstruc- 
tion of the mass and the cross-section assuming an annihilation into bb. The true point is indicated 
by a symbol; the contours delimit the regions that Ams-02 would identify, at the indicated C.L. and 
given the indicated assumptions on the background. Middle Panel: Reconstruction of the mass and 
the cross-section without knowing a priori the annihilation channel. The true point is indicated by the 
symbol; the other points and the contours indicate the best fits and the 90% C.L. regions that Ams-02 
would identify for other annihilation channels. Right Panels: Reconstruction of the mass and the 
branching ratio for a fixed cross-section and annihilation into bb and W + W~ . The symbols represent 
the model used to produce the mock data. 
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Fixed annihilation channel 


mass niDM 


cross-section (av) 


Ax 2 with respect 


(true signal: 100% bb) 


[GeV] 


[cm 3 s _1 ] 


to a pure background 


XX -> bb 


84.9 


2.7- 10~ 25 


-10.5 


xx ->• w+w- 


84.8 


6.6 • 10~ 26 


-10.3 


XX ^ 77 


57.4 


1.0 • 10~ 24 


-5.7 


Fixed cross-section (av) [cm 3 s _1 ] 


mass ?TirjM 


relative 


Ax' 2 with respect 


(true signal: 70% bb + 30% W + W~) 


[GeV] 


branching ratio 


to a pure background 


2.7 x 10~ 26 


84.9 


0.7 


-5.4 



Table 2: Best fit reconstructions of a signal from a 85 GeV DM candidate with (av) = 2.7 x 
10~ 26 cm 3 s _1 (model k). See text for details. 



Fixed annihilation channel 


mass raoM 


cross-section (av) 


Ax' 2 with respect 


(true signal: 100% bb) 


[GeV] 


[cm 3 s _1 ] 


to a pure background 


XX -> bb 


300 


1.5 • 10" 25 


-21.0 


xx -> w+w- 


240 


1.9 • 10~ 25 


-19.7 


XX ->• 77 


169 


4.8 • 10' 24 


-9.8 


->• 


447 


7.6 • 10~ 23 


-19.2 


Fixed cross-section (av) [cm 3 s~ 1 ] 


mass moM 


relative 


Ax 2 with respect 


(true signal: 40% bb + 60% W + W~ ) 


[GeV] 


branching ratio 


to a pure background 


1.5 x 10~ 25 


300 


0.4 


-13.0 



Table 3: Best fit reconstructions of a signal from a 300 GeV DM candidate with (av) = 1.5 x 
10~ 25 cm 3 s" 1 (model -k). See text for details. 



Fixed annihilation channel 


mass niDM 


cross-section (av) 


Ax 2 with respect 


(true signal: 100% bb) 


[GeV] 


[cm 3 s _1 ] 


to a pure background 


XX -> bb 


999 


5 • 10" 25 


-15.7 


xx ->• w+w- 


886 


5.8 • 10" 24 


-15.1 


XX -> 77 


765 


2.5 -lO" 22 


-14.5 


xx ->• 


1711 


1.1 • io- 22 


-14.8 


Fixed cross-section (av) [cm 3 s _1 ] 


mass riiDM 


relative 


Ax' 2 with respect 


(true signal: 60% bb + 40% W+W~ ) 


[GeV] 


branching ratio 


to a pure background 


5 x lO" 25 


999 


0.6 


-15.3 



Table 4: Best fit reconstructions of a signal from a 1 TeV DM candidate with (av) = 5x 10 25 cm 3 s 
(model m). See text for details. 
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hypothesis. 2 We see that, while the bb channel does still obtain the overall best fit, all the other 
channels also fit the data similarly well, so that it is not practically possible to reconstruct the 
'true' configuration. In addition, we remind that a channel like tt (when kinematically open) 
is practically indistinguishable from bb and therefore constitutes an additional degeneracy. 

The Ax 2 numbers allow to draw another conclusion: all channels for the * and ■ models 
(with the possible exception of the borderline case of the 77 channel in the ▲ model) improve 
sufficiently the x 2 that an excess can quite clearly be identified. This might be surprising if, 
for instance, judging by eye from the right panel of fig. 1. However, as long as a power-law 
assumption (with the assumed uncertainty band) for the background is justified, this statistical 
conclusion is sound. If one allows for the possibility that the background shows features mim- 
icking the hump produced by DM (or that it might consist of a superposition of power-laws, 
or that the uncertainty band in much larger...), then of course these conclusions evaporate, as 
warned in Sec. 2.2. On the other hand, for the ▲ model the improvement of x 2 over the pure 
background becomes quite thin: an actual reconstruction will be very challenging. 

Branching Ratio 

Finally, we produce mock data with a mixed annihilation into bb and into W + W~ and we 
investigate whether it is possible to reconstruct the Branching Ratio. The results in the right 
panels of fig. 5 show that it is not possible. While the determined best fit falls on top or very 
close to the true value (see the last lines of the tables 2, 3, 4), the 90% and 99% C.L. regions 
span the entire range of BR. This is due to the similarity of the bb and W + W~ spectra and it 
reinforces the conclusions of the analyses above. In particular, these conclusions highlight that 
it will be important to carry on a multi-messenger analysis (i.e. including signals in 7 rays, 
neutrinos etc) in order to break the degeneracies and obtain a convincing reconstruction. 

4 Conclusions 

In this paper we have explored the potential of the current and future measurements of the 
flux of CR antiprotons to constrain Dark Matter. We have analyzed, in a model independent 
way, several annihilation and decay channels, a large range of masses and assessed the impact 
of astrophysical uncertainties. Our main results can be schematized as follows. 

— The current p constraints using PAMELA measurements are very strong both for annihi- 
lating (see fig. 2 left) and decaying (see fig. 3 left) Dark Matter, for the hadronic chan- 
nels. Adopting fiducial choices for the halo profile and propagation parametres ('Einasto, 
MED'), they are as strong as (or even slightly stronger than) the most stringent gamma 
ray constraints from the FERMI satellite. 

— The astrophysical uncertainty band associated to different choices for the halo profile and 
the propagation parameters, however, spans between one and two orders of magnitude 
(see fig. 2 right and fig. 3 right). 

2 More precisely, what we do is the following. We generate the mock data using a 'fixed background + signal' 
flux. We then try and interpret them in terms of a 'variating background' with free A and p: the best of such 
variating backgrounds has a \ 2 of Xo- We then add back the signal, assuming different annihilation channels, 
and we see how much the \ 2 improves with respect to \o- These improvements are the numbers we quote in 
the last column of the tables. 
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— The upcoming Ams-02 measurements (assuming 1 year of accumulated data and a real- 
istic acceptance for the experiment) have the power to improve the constraints, assuming 
that no signal is seen, by slightly less than one order of magnitude (see fig. 4). In partic- 
ular, Ams-02 will probe the thermal value of the annihilation cross section in the range 
todm — 30 — 200 GeV for the bb (or tt) channel. 

— If a convincing excess is measured by Ams-02 and it is attributed to annihilating Dark 
Matter, the data will allow to somewhat constrain the underlying DM properties, but 
a full reconstruction will be quite challenging (see fig. 5). The main reason is that the 
p spectra from different channels are quite similar one to another (as a result of the 
hadronization process smoothing out differences) and not very different from the shape 
of the pure background, so that degeneracies are impossible to break. 

— Restricting to the 'traditional' hadronic annihilation channels (bb, ti, W + W~ , ZZ), 
in the most favorable case in which moM ~ few hundreds GeV (and therefore the 
signal falls in the energy range best measured by Ams-02) it will be possible to 
reconstruct the DM mass within 50%. For smaller values of tbdm (below ~ 100 GeV) 
or larger ones (above ~ 1 TeV) the reconstruction capability worsens somewhat. 

— Still restricting to the hadronic channels, the determination of the annihilation cross 
section will be possible within an order of magnitude independently of the mass 
range. 

— If other annihilation channels are allowed, the data will not allow to discriminate 
among them and the determination of the annihilation cross section becomes impos- 
sible. 

— If several of the hadronic channels are open at the same time, it will not be possible 
to determine the relative branching ratios. 

In summary, we find that antiprotons are a very relevant tool to constrain Dark Matter annihi- 
lation and decay, on a par with gamma rays for the hadronic channels. Current Pamela data 
and especially upcoming Ams-02 data allow to probe large regions of the parameter space. On 
the other hand, using p to reconstruct the DM properties in case of positive detection will be 
very challenging, except for favorable scenarios. 
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